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Abstract. The behavior of a polyelectrolyte adsorbed on a charged substrate of high-dielectric constant is studied by both 
Monte-Carlo simulation and analytical methods. It is found that in a low enough ionic strength medium, the adsorption 
transition is first-order where the substrate surface charge still keeps repulsive. The monomer density at the adsorbed surface 
is identified as the order parameter. It follows a linear relation with substrate surface charge density because of the electrostatic 
boundary condition at the charged surface. During the transition, the adsorption layer thickness remains finite. A new scaling 
law for the layer thickness is derived and verified by simulation. 



INTRODUCTION 

The problem of polymer adsorption on an attractive sur- 
face has drawn considerable interest due to its relation 
to surface effects in critical phenomena and practical im- 
portance in technology and biology. It is well established 
that the adsorption transition is continuous if its attrac- 
tion on the surface is short-rang ed[l]. On the other hand, 
long-ranged electrostatic interactions in polyelectrolyte 
systems pose many challenging theoretical problems. 
Recently the macroion adsorption on an electrostatically 
attractive interface |2j |3[ and the associated charge in- 
version phenomena of adsorbed polyelectrolytes @, |5| 
acquire lots of attention. 

Previous analytical approaches impose the continuity 
of monomer density across the charged surface, the sur- 
face monomer density is then set to zero |6, 7]. The poly- 
electrolyte adsorption problem becomes the competition 
between the electrostatic force from charged surface and 
the configurational entropy of the polyelectrolyte itself. 
The transition has been shown to be continuous. How- 
ever, we found that the above treatment is not quite ad- 
equate since the electrostatic boundary condition at the 
charged surface is not faithfully respected. Furthermore, 
the electrostatic boundary condition may induce different 
physics in the presence of a substrate under the charged 
surface. In this paper, we study the adsorption of a poly- 
electrolyte on a high-dielectric substrate at low ionic 
strength (e.g. aqueous solution with a metal substrate) 
in which image charge interaction is attractive. The ad- 
sorption transition occurs where the surface charges are 
repulsive instead of the attractive case that were usually 
studied. The problem is tackled by Monte-Carlo simula- 
tions and analytical methods taking full account of ap- 
propriate boundary conditions. It is found that the order 



of adsorption transition, the physical mechanism, and the 
scaling behavior are all different from those of the attrac- 
tive surfaces. 



MODEL 

A polyelectrolyte carrying positive charges is immersed 
in a medium (z > 0) of dielectric constant e. At z = 
there is an impenetrable surface of uniformly surface 
charge density a. Below that (z < 0), it is a substrate 
of dielectric constant e'. Denote the charge on a polymer 
segment ds by qods, the Hamiltonian is written as 



3k B T 
2l 2 



o V ds 



1 rN rN 

- ds ds'{ 

2 Jo Jo 



WFW)\ 



- T' (2 - C\y ) — — —— } 

\r(s)-7(s')\ 



dsK-^-* 7 ^ 







rN rN 

-CO I ds ds'8(r(s)-r(s')) 
Jo Jo 



(1) 



where s is the variable to parametrize the chain, 
Iq the bare persistence length, and K~ l the De- 
bye screening length. r(s) = (x(s),y(s),z(s)), 
r'(i') = (x(s'),y(s'),—z(s')) are the positions of 
the monomers and their electrostatic images, re- 
spectively, r = ql/e, r' = r(e' - e)/(e' + e), and 
h = 4-nqoG / (e' + e) are the coupling parameters gov- 
erning the strengths of Coulomb interactions among 
the monomers themselves, between the polymer and 
its image, and between the polymer and the charged 



surface, respectively. In good solvent regime, co > in 
the last term. 



MONTE-CARLO RESULTS 

The continuum model is discretized to perform Monte- 
Carlo simulations. r(s) is replaced by a chain of beads F, 
(/ = l,...,N) with hard-core excluded volume of finite 
radius a. Polymer lengths up to = 120 are employed. 
The units of length and energy are 2a and q\/2ea, re- 
spectively. The ratios e'/e are studied from 2 to 12.5 
(aqueous solution with a metallic substrate). Runs up to 
10 9 MC steps are performed. 

The adsorption behavior can be characterized by nor- 
malized monomer density p(z). p a = p(z = a), the prob- 
ability density of monomers adsorbed on the substrate, is 
chosen as the order parameter to describe the adsorption 
transition. p u as a function of the surface charge density 
a for various e'/e > 1 is shown in Fig[2 It is seen that 
p a vanishes abruptly when a increases up to its thresh- 
old value (7 t indicating a first-order transition. At low 
enough ionic strength that K ' is much larger than ab- 
sorption layer thickness, a t > for e'/ e > 1. The discon- 
tinuous drop of p a across the transition decreases to zero 
as e'/e — * 1. We have also verified that the sharp jump in 
energy (latent heat) across the transition is proportional 
to N as expected for a first-order transition. Same results 
are obtained for larger fc -1 . 

Furthermore, p a is linear in a with the slope depend- 
ing on e'/e. Such a linear relation between p a and a can 
be understood from the electrostatic boundary conditions 
that the system has to satisfy. The electric potential <p (z) 
in the neighborhood of z = boundary obeys 
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where a p is the polarization surface charge density in- 
duced by the polyelectrolyte only. Notice that c p depends 
only on e'/e. It is independent of a in the adsorbed 
regime near the transition. The reason is based on the 
electric blob picture which will be explained in next sec- 
tion. If one treats the polyelectrolyte as a marcomolecule 
with a well-defined surface, its surface charge density at 
z = a should be proportional to the monomer density p„. 
It also applies to the electric field in the z < region. One 
have 
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where K > is the proportionality constant. Applying 
the electric field continuity from z — + to z = cC , and 




FIGURE 1. Monte-Carlo results for the normalized 
monomer density at the surface, p fl , as a function of surface 
charge density a (in unit of qa/Aa 2 ) for different e'/e at 
K =25. The fitted straight lines are terminated at their 
adsorption transition points. The vertical dashed lines are 
drawn as guides to the eyes. Inset: The polarization surface 
charge density due to the charged polymer, | o" p | , as a function 
of e'/e. Op saturates as e'/e — » °°. The sign of <7 p is opposite 
to go and is thus negative. The solid curve is fitted from Eq.{6) 
with a poly = 0.118. 



using Eqs.(|2|i-@, one get the linear behavior 
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Notice that K and <7 p depends only on e'/e. As shown 
in FigQ the slope K decreases monotonically with e'/e. 
Lines with different values of e'/ e intersect at a common 
point implies the naive intuition that a more polarizable 
substrate gives a stronger adsorption is not always true. 

Without substrate surface charge, set a = in Eq.0, 
one get 
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The dependence on e'/e can be considered as a flat sur- 
face of surface charge density <7 po i y = p a \a=o/K located 
above the substrate. K are obtained by linear fittings in 
FigH and hence (7 p as a function of e'/e is plotted in 
the inset. (7 p increases from zero to a saturated value as 
e'/ e — > oo. Presumably, Eq.© suggests that it is possible 
to experimentally detect the sign and magnitude of the 
charge of the polyelectrolyte by measuring the substrate 
polarization. 



ANALYTICAL RESULTS 

At low ionic strength, a polyelectrolyte can be treated 
as electric blobs arranged longitudinally. The rod-like 



polyelectrolyte tends to lie down on the substrate so 
as to lower its energy. It is thus rigid (rod-like) in x-y 
plane, but still flexible in the z-direction. In the adsorbed 
regime, decreasing a (chain is attracted more by the sur- 
face) would cause the rearrangement of the electric blobs 
such that the fluctuations in the z-direction is reduced 
(and hence the layer thickness decreases), but the fluctu- 
ations in the other two directions and the size of the blobs 
are basically unchanged. The effective in-plane surface 
charge distribution of the polyelectrolyte, and hence <7 p , 
is not affected by a. The excluded volume effect is ig- 
nored because it takes almost no effect in z-direction. 
The effect from self-electrostatic interaction can be ab- 
sorbed by renormalising the bare persistence length Iq to 

Because the monomer would feel the strongest attrac- 
tion from its direct image around the adsorption regime, 
the r'-term in Eq.Q is approximated by the interaction 
of every monomer and its corresponding image only. The 
residual attraction from the images of other monomers 
could be absorbed into the coupling parameter F' by 
renormalising qo to q. The partition function is reduced 
to 
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r(s) is transformed to p(r) = ^ J ds8(r — r{s)) by in- 
troducing an auxiliary field. Ground state dominance in 
large-Af limit is then applied. By variational principle, 
one obtains the Edwards-Schrodinger equation, 
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where Sq acts as a Lagrange multiplier to enforce the 
constraint of the ground state wavefunction normaliza- 
tion. The monomer density is given by p(z) = |y/(z)| 2 . 
Eq.© also describes a quantum particle at its ground 
state moving under a combined potential of a Id screened 
Coulomb attraction and an almost linear potential. How- 
ever, the boundary condition expressed by Eq.Q is dif- 
ferent from that of the hard- wall one i//"| s — usually em- 
ployed for the quantum particle. Instead y/| s = yfp2 ^ 
for the present problem implies that the steric force felt 
by the polyelectrolyte from the charged surface should be 
modified (§]. Setting y/| s = in previous studies 

is not completely correct. 

During the adsorption, the rod-like polyelectrolyte 
tends to lie down on the charged surface. The thickness 
of the adsorption layer is of the same order of the gy- 
ration radius in z-direction. At low ionic strength corre- 
sponding to Debye length much larger than layer thick- 
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FIGURE 2. Effective Potential felt by the polyelectrolyte 



ness, the polyelectrolyte can only feel the potential bar- 
rier formed by V(z) as shown in Fi^2] Analytically, the 
original potential V (z) in Eq.(|8} is replaced by 
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where z b i is chosen such that V'(zbr) = 0, and V(zbr) is 
the barrier height. In the limit of a = K = 0, the analytic 
solution exists in the form 
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where ^ / 2 is the Whittaker's notation of the confluent 



hypergeometric function 1 10], and X is the least value 
satisfying the boundary condition. Bound state exists for 
arbitrary e'/e > 1. It implies Oi > at low enough ionic 
strength. 

No exact solution exists for a, K > in general but 
one can analyze it around the transition. In the region 
of interest, the Coulomb term dominates over the almost 
linear term. One can approximate the binding energy 
by Eq.dim. The polyelectrolyte undergoes a de-sorption 
transition when the binding energy is equal to barrier 
height V(zbr)- After some algebra, it shows a t ~ (e'/e — 
1) for e'/e > 1 and a t - (e'/e - l) 3 for e'/e - 1. This 
analytic result is consistent with simulation as shown in 
FigEt. 

An approximate solution for density profile p (z) for 
a > can be studied by the variational wavefunction 

W(z) = V7^(l +Ai«(z-a))e-2<- fl ) (11) 

where a -1 is the decay length, is positive because the 
wavefunction is restricted to be nodeless. a and jj. are not 
independent but related via the wavefunction normaliza- 
tion condition. The inverse decay length is calculated to 
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FIGURE 3. (a) Simulation results for the surface charge 
density at the transition, a t , as a function of e'/e in logarithmic 
scale at ic -1 = 25. The straight lines indicate slopes of 1 and 
3 as suggested in the text, (b) Simulation results for the decay 
length a -1 (which is proportional to the layer thickness) as a 
function of a (in unit of qo/4a 2 ) for e'/e = 12.5, K" -1 = 25. 
The straight line is a linear fit to points near the transition. 
o"t = 0. 102 in this case. a~ 1 is obtained from exponential fitting 
to the tail of corresponding density profile. 

be 

3BT' 
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where the leading term is independent of <T. Near the 
transition, the decay length a -1 and hence the adsorp- 
tion layer thickness D ~ a -1 remains finite. Combining 
Eqs.© and 1121 . one get the scaling behavior 

D t -Z)~(af 1 -a _1 )~(a t -a) (13) 

where D t is the threshold layer thickness, a -1 as a func- 
tion of a from simulations is shown in Fig0). 

It would be beneficial to compare with the case that 
a polyelectrolyte is adsorbed onto an attractive charged 
substrate of low-dielectric constant (e.g. DNA in aque- 
ous solution onto a charged lipid membrane). Its asymp- 
totic solution for z — ► +°° to Eq.(|8} reproduces the usual 
scaling/)^ a -1 ~ [4]. The adsorption onto a low- 

dielectric substrate is continuous. The order parameter p„ 
vanishes continuously across the transition, and the layer 
thickness swells to infinity as the polyelectrolyte is de- 
sorbed 111 ill . 



is the Bjerrum length). Similarly, the strongly charged 
surface ofbare charge density larger than K/(nls) is just 
renormalized back to k/(jcIb) UM - 

Although we focus on the adsorption of a single poly- 
electrolyte, our results may provide a starting point to 
study the charge inversion and multi-layer adsorption. 
(3- At low ionic strength, polyelectrolytes are adsorbed 
in a multi-layer structure because of strong Coulomb 
repulsion. Each layer is composed of the parallel Id 
Wigner crystal 1 15]. From our physical picture of a single 
polyelectrolyte adsorption, the upper bound for the thick- 
ness of the multi-layer structure is Zbr ~ o~ 1 ' 2 {s l /e — 
l) 1 / 2 . This suggests one can easily control just one sin- 
gle layer adsorbed onto high-dielectric substrate by tun- 
ing the surface charge density. Rigorous treatment based 
on this physical picture will be elaborated elsewhere. 
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DISCUSSION 

A strongly charged polyelectrolyte immersed in a salt 
solution will attract oppositely charged ions to condense 
until its effective charge density reaches the Manning 
threshold 11211 . This means that one can just renormalize 
qo in our system to lea / 1^ if qo is larger than lea /Zg (Is 



